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1 INTRODUCTION 

Astrophysical and cosmological arguments and observations have become part 
of the main-stream methodology to obtain empirical information on existing or 
hypothetical elementary particles and their interactions. The "heavenly laborato- 
ries" are complementary to accelerator and non-accelerator experiments, notably 
at the "low-energy frontier" of particle physics, which includes the physics of neu- 
trinos and other weakly interacting low-mass particles such as the hypothetical 
axions, novel long-range forces, and so forth. 

The present review is dedicated to stars as particle-physics laboratories, or 
more precisely, to what can be learned about weakly interacting low-mass par- 
ticles from the observed properties of stars. The prime argument is that a hot 
and dense stellar plasma emits low-mass weakly interacting particles in great 
abundance. They subsequently escape from the stellar interior directly, without 
further interactions, and thus provide a local energy sink for the stellar medium. 
The astronomically observable impact of this phenomenon provides some of the 
most powerful limits on the properties of neutrinos, axions, and the like. 

Once the particles have escaped they can decay on their long way to Earth, 
allowing one to derive interesting limits on radiative decay channels from the 
absence of unexpected x- or 7-ray fluxes from the Sun or other stars. 

Finally, the weakly interacting particles can be directly detected at Earth, thus 
far only the neutrinos from the Sun and supernova (SN) 1987A, allowing one to 
extract important information on their properties. 

The material covered here has been reviewed in 1990, with a focus on axion 
limits, by Turner [ |TJ and by Raffelt [ [|, and a very brief "Mini-Review" was 
included in the 1998 edition of the Review of Particle Physics [ [3]. My 1996 
book Stars as Laboratories for Fundamental Physics [ |I| treats these topics in 
much greater detail than is possible here. The present chapter is intended as a 
compact, up-to-date, and easily accessible source for the most important results 
and methods. 

The subject of "Particle Physics from Stars" is broader than both my expertise 
and the space available here. I will not touch on the solar neutrino problem and 
its oscillation interpretation — this is a topic unto itself and has been extensively 
reviewed by other authors, for example [ [5|, ||, [?]]. 

The high densities encountered in neutron stars make them ideal for studies and 
speculations concerning novel phases of nuclear matter (e.g. meson condensates 
or quark matter), an area covered by two recent books [ ||, [J. Quark stars are 
also the subject of an older review [ |l(]] and are covered in the proceedings of two 
topical conferences [ 11, ^] 



Certain grand unified theories predict the existence of primordial magnetic 
monopoles. They would get trapped in stars and then catalyze the decay of 
nucleons by the Rubakov-Callan effect. The ensuing anomalous energy release 
is constrained by the properties of stars, in particular neutron stars and white 
dwarfs, a topic that has been reviewed a long time ago [ [l3|]. It was re-examined, 
and the limits improved, in the wake of the discovery of the faintest white dwarf 
ever detected which puts restrictive limits on an anomalous internal heat source [ 
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Weakly interacting massive particles (WIMPs), notably in the guise of the 
supersymmetric neutralinos, are prime candidates for the cosmic dark matter. 
Some of them would get trapped in stars, annihilate with each other, and produce 
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a secondary flux of high-energy neutrinos. The search for such fluxes from the 
Sun and the center of the Earth by present-day and future neutrino telescopes is 
the "indirect method" to detect galactic particle dark matter, an approach which 



is competitive with direct laboratory searches — see [ 15] for a review. 

Returning to the topics which are covered here, Sections |^-|] are devoted to a 
discussion of the main stellar objects that have been used to constrain low-mass 
particles, viz. the Sun, globular-cluster stars, compact stars, and SN 1987A. In 
Sections |5|-|7| the main constraints on neutrinos, axions, and novel long-range 



forces are summarized. Section |8| is given over to brief concluding remarks. 



2 THE SUN 



2.1 Basic Energy-Loss Argument 

The Sun is the best-known star and thus a natural starting point for our survey 
of astrophysical particle laboratories. It is powered by hydrogen burning which 
amounts to the net reaction 4p + 2e~ — > 4 He + 2i^ e + 26.73 MeV, giving rise to 
a measured v e flux which now provides one of the most convincing indications 
for neutrino oscillations [ g, g, 0. Instead of neutrinos from nuclear processes 
we focus here on particle fluxes which are produced in thermal plasma reactions. 
The photo neutrino process 7 + e~ — ► e~ + vv is a case in point, as is the 
production of gravitons from electron bremsstrahlung. The solar energy loss 
from such standard processes is small, but it may be large for new particles. To 
be specific we consider axions (Sec. |6|) which arise in a variety of reactions, and in 
particular by the Primakoff process in which thermal photons mutate into axions 
in the electric field of the medium's charged particles (Fig. [T]). In Sec. 6.2.1 we 
will discuss direct search experiments for solar axions, while here we focus on 
what is the main topic of this review, the backreaction of a new energy loss on 
stars. 



Figure 1: Primakoff production of axions in the Sun. 



The Sun is a normal star which supports itself against gravity by thermal 
pressure, as opposed to degenerate stars like white dwarfs which are supported 
by electron degeneracy pressure. If one pictures the Sun as a self-gravitating 
monatomic gas in hydrostatic equilibrium, the "atoms" obey the virial theorem 
(-^kin) = —\ {-f'grav}- The most important consequence of this relationship is that 
extracting energy from such a system, i.e. reducing the total energy (i?km) + 
(.Egrav), leads to contraction and to an increase of (-Ekin)- Therefore, all else 
being equal, axion losses lead to contraction and heating. The nuclear energy 
generation rate scales with a high power of the temperature. Therefore, the 
heating implied by the new energy loss causes increased nuclear burning — the 
star finds a new equilibrium configuration where the new losses are compensated 
by an increased rate of energy generation. 

The main lesson is that the new energy loss does not "cool" the star; it leads to 
heating and an increased consumption of nuclear fuel. The Sun, where energy is 
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transported from the central nuclear furnace by radiation, actually overcompen- 
sates the losses and brightens, while it would dim if the energy transfer were by 
convection. Either behavior is understood by a powerful "homology argument" 
where the nonlinear interplay of the equations of stellar structure is represented 
in a simple analytic fashion [ [H]] . 

The solar luminosity is well measured, yet this brightening effect is not observ- 
able because all else need not be equal. The present-day luminosity of the Sun 
depends on its unknown initial helium mass fraction Y; in a solar model Y has to 
be adjusted such that Lq = 3.85 x 10 33 erg s~ x is reproduced after 4.6 x 10 9 years 
of nuclear burning. For solar models with axion losses the required presolar he- 
lium abundance Y as a function of the axion-photon coupling constant g a ^ is 
shown in Table |l]. The axion luminosity L a is also given as well as the central 
helium abundance Y c , density p c , and temperature T c of the present-day Sun. 

Even axion losses as large as Lq can be accommodated by reducing the preso- 
lar helium mass fraction from about 27% to something like 23% [ |17|, [lql . The 
"standard Sun" has completed about half of its hydrogen-burning phase. There- 
fore, the anomalous energy losses cannot exceed approximately Lq or else the 
Sun could not have reached its observed age. Indeed, for gio = 30 no consistent 



present-day Sun could be constructed for any value of Y { 18 1. The emission rate 
of other hypothetical particles would have a different temperature and density 
dependence than the Primakoff process, yet the general conclusion remains the 
same that a novel energy loss must not exceed approximately Lq. 



2.2 Solar Neutrino Measurements 

This crude limit is improved by the solar neutrino flux which has been measured 
in five different observatories with three different spectral response characteristics, 
i.e. by the absorption on chlorine, gallium, and by the water Cherenkov technique. 
The axionic solar models produce larger neutrino fluxes; in Table [| we show the 
expected detection rates for the CI, Ga, and H2O experiments relative to the 
standard case. For g\$ < 10 one can still find oscillation solutions to the observed 
v e deficit, but larger energy-loss rates appear to be excluded [ |i~7f . 

Once the neutrino oscillation hypothesis has been more firmly established and 
the mixing parameters are better known, the neutrino measurements may be used 
to pin down the central solar temperature, allowing one to constrain novel energy 
losses with greater precision. For now it appears safe to conclude that the Sun 



Table 1: Solar-model parameters and relative detection rates in the CI, Ga, and 
water neutrino observatories as a function of the axion-photon coupling constant 
5 a7 /(10~ 10 GeV" 1 ) according to Ref. 
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does not emit more than a few tenths of L Q in new forms of radiation 



2.3 Helioseismology 

Over the past few years the precision measurements of the solar p-mode frequen- 
cies have provided a more reliable way to study the solar interior. For example, 
the convective surface layer is found to reach down to 0.710-0.716 Rq [ [L9|, the 
helium content of these layers to exceed 0.238 [ 2C]. Gravitational settling has 



reduced the surface helium abundance by about 0.03 so that the presolar value 
must have been at least 0.268, in good agreement with standard solar models. 
The reduced helium content required of the axionic solar models in Table |l| dis- 
agrees significantly with this lower limit for gio > 10. 

One may also invert the p-mode measurements to construct a "seismic model" 
of the solar sound-speed profile, e.g. [ pi}] . Ah 1 modern standard solar models 
agree well with the seismic model within its uncertainties (shaded band in Fig. |2|) 
which mostly derive from the inversion method itself, not the measurements. The 
difference between the sound-speed profile of a standard solar model and those 
including axion losses are also shown in Fig. |2} For gio > 10 the difference is 
larger than the uncertainties of the seismic model, implying a limit 

g ay < 10 x 10~ 10 GeV" 1 . (1) 

Other cases may be different in detail, but it appears safe to assume that any 
new energy- loss channel must not exceed something like 10% of Lq. 




-0.02 r 



-0.03 r V. 

-0.04 F , 
0.0 



0.2 



0.4 



0.6 



0.8 



r/Ro 



Figure 2: Fractional difference in sound-speed profiles of solar models with axion 
losses compared to the reference model [ [l?]] . The shaded area is the uncertainty 
of the seismic model [pOfl. The axion-photon coupling constant was <7io=4.5 (solid 
line), 10 (short-dashed), 15 (dash-dotted), 20 (dash-dot-dot-dotted). 



2.4 "Strongly" Interacting Particles 

Thus far we have assumed that the new particles couple so weakly that they 
escape from the stellar interior without further interactions, in analogy to neutri- 
nos or gravitons. They emerge from the entire stellar volume, i.e. their emission 
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amounts to a local energy sink for the stellar plasma. But what if the particles 
interact so strongly that their mean free path is less than the solar radius? 

The impact of such particles on a star compares to that of photons, which 
are also "trapped" by their "strong" interaction. Their continuous thermal pro- 
duction and re-absorption amounts to the net transfer of energy from regions 
of higher temperature to cooler ones. In the Sun this radiative form of energy 
transfer is more important than conduction by electrons or convection, except in 
the outer layers. A particle which interacts more weakly than photons is more 
effective because it travels a larger distance before re-absorption — the ability to 
transfer energy is proportional to the mean free path. The properties of the Sun 
roughly confirm the standard photon opacities, so that a new particle would have 



to interact more strongly than photons to be allowed [ 22, |2 



Therefore, contrary to what is sometimes stated in the literature, a new particle 
is by no means allowed just because its mean free path is less than the stellar 
dimensions. The impact of a new particle is maximal when its mean free path 
is of order the stellar radius. Of course, usually one is interested in very weakly 
interacting particles so that this point is moot. 



3 LIMITS ON STELLAR ENERGY LOSSES 



3.1 Globular- Cluster Stars 
3.1.1 Evolution of Low-Mass Stars 

The discussion in the previous section suggests that the emission of new weakly 
interacting particles from stars primarily modifies the time scale of evolution. 
For the Sun this effect is less useful to constrain particle emission than, say, the 
modified p-mode frequencies or the direct measurement of the neutrino fluxes. 
However, the observed properties of other stars provide far more restrictive limits 
on their evolutionary time scales so that anomalous modes of energy loss can be 
far more tightly constrained. We begin with globular-cluster stars which, together 
with SN 1987A, are the most successful example of astronomical observations that 
provide nontrivial limits on the properties of elementary particles. 




Figure 3: Globular cluster M3. (Image courtesy of Palomar/Caltech.) 
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Figure 4: Color-magnitude diagram for the globular cluster M3, based on the 
photometric data of 10,637 stars [ Vertically is the brightness in the visual 
(V) band, horizontally the difference between B (blue) and V brightness, i.e. a 
measure of the color and thus surface temperature, where blue (hot) stars lie 
toward the left. The classification for the evolutionary phases is as follows [ 
|25fl . MS (main sequence): core hydrogen burning. BS (blue stragglers). TO 
(main-sequence turnoff): central hydrogen is exhausted. SGB (subgiant branch): 
hydrogen burning in a thick shell. RGB (red-giant branch): hydrogen burning 
in a thin shell with a growing core until helium ignites. HB (horizontal branch): 
helium burning in the core and hydrogen burning in a shell. AGB (asymptotic 
giant branch): helium and hydrogen shell burning. P-AGB (post-asymptotic 
giant branch): final evolution from the AGB to the white-dwarf stage. 

Our galaxy has about 150 globular clusters such as M3 (Fig. |||) which are gravi- 
tationally bound systems of up to a million stars. In Fig. [| the stars of the cluster 
M3 are arranged according to their color or surface temperature (horizontal axis) 
and brightness (vertical axis) in the usual way, leading to a characteristic pattern 
which allows for rather precise tests of the theory of stellar evolution, and notably 
for quantitative measurements of certain evolutionary time scales. Globular clus- 
ters are the oldest objects in the galaxy and thus almost as old as the universe. 
The stars in a given cluster all formed at about the same time with essentially the 
same chemical composition, differing primarily in their mass. Because more mas- 
sive stars evolve faster, present-day globular-cluster stars are somewhat belowQ 
1M.Q so that we are concerned with low-mass stars {M. < 2M.q). Textbook 
expositions of stellar structure and evolution are [ |26|, p7[ . 

Stars begin their life on the main sequence (MS) where they burn hydrogen 
in their center. Different locations on the MS in a color-magnitude diagram 
like Fig. |] correspond to different masses, with more massive stars shining more 

lr The letter M denotes stellar masses with IMq = 2 x 10 33 g the solar mass. The letter M 
is traditionally reserved for the absolute stellar brightness (in magnitudes or mag). The total 
or bolometric brightness is defined as A/boi = 4.74 — 2.5 log 10 (L / Lq) , with the solar luminosity 
L Q = 3.85 x 10 33 erg s" 1 . 
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brightly. When central hydrogen is exhausted the star develops a degenerate 
helium core, with hydrogen burning in a shell. Curiously, the stellar envelope 
expands, leading to a large surface area and thus a low surface temperature (red 
color) — they become "red giants." The luminosity is governed by the gravita- 
tional potential at the edge of the growing helium core so that these stars become 
ever brighter: they ascend the red-giant branch (RGB). The higher a star on the 
RGB, the more massive and compact its helium core. 

The core grows until about 0.5 A4® when it has become dense and hot enough to 
ignite helium. The ensuing core expansion reduces the gravitational potential at 
its edge and thus lowers the energy production rate in the hydrogen shell source, 
dimming these stars. Helium ignites at a fixed core mass, but the envelope mass 
differs due to varying rates of mass loss on the RGB, leading to different surface 
areas and thus surface temperatures. These stars thus occupy the horizontal 
branch (HB) in the color-magnitude diagram. In Fig. || the HB turns down on 
the left (blue color) where much of the luminosity falls outside the V filter; in 
terms of the total or "bolometric" brightness the HB is truly horizontal. 

Finally, when helium is exhausted, a degenerate carbon-oxygen core develops, 
leading to a second ascent on what is called the asymptotic giant branch (AGB). 
These low-mass stars cannot ignite their carbon-oxygen core — they become white 
dwarfs after shedding most of their envelope. 

The advanced evolutionary phases are fast compared with the MS duration 
which is about 10 10 yr for stars somewhat below IMq. For example, the ascent 
on the upper RGB and the HB phase each take around 10 s yr. Therefore, the 
distributions of stars along the RGB and beyond can be taken as an "isochrone" 
for the evolution of a single star, i.e. a time-series of snapshots for the evolution of 
a single star with a fixed initial mass. Put another way, the number distribution 
of stars along the different branches are a direct measure for the duration of the 
advanced evolutionary phases. The distribution along the MS is different in that 
it measures the distribution of initial masses. 

3.1.2 Core Mass at Helium Ignition 

Anomalous energy losses modify this picture in measurable ways. We first con- 
sider an energy-loss mechanism which is more effective in the degenerate core of a 
red giant before helium ignition than on the HB so that the post-RGB evolution 
is standard. Since an RGB-star's helium core is supported by degeneracy pres- 
sure there is no feedback between energy-loss and pressure: the core is actually 
cooled. Helium burning (3 4 He — > 12 C) depends very sensitively on temperature 
and density so that the cooling delays the ignition of helium, leading to a larger 
core mass A4 C , with several observable consequences. 

First, the brightness of a red giant depends on its core mass so that the RGB 
would extend to larger luminosities, causing an increased brightness difference 
AM^g between the HB and the RGB tip. Second, an increased A4 C implies an 
increased helium-burning core on the HB. For a certain range of colors these stars 
are pulsationally unstable and are then called RR Lyrae stars. From the measured 
RR Lyrae luminosity and pulsation period one can infer M c on the basis of their 
so-called mass-to-light ratio A. Third, the increased A4 C increases the luminosity 
of RR Lyrae stars so that absolute determinations of their brightness Mrr allow 
one to constrain the range of possible core masses. Fourth, the number ratio R 
of HB stars vs. RGB stars brighter than the HB is modified. 
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These observables also depend on the measured cluster metallicity as well as the 
unknown helium content which is usually expressed in terms of Y env , the envelope 
helium mass fraction. Since globular clusters formed shortly after the big bang, 
their initial helium content must be close to the primordial value of 22-25%. 
Y env should be close to this number because the initial mass fraction is somewhat 
depleted by gravitational settling, and somewhat increased by convective dredge- 
up of processed, helium-rich material from the inner parts of the star. 

An estimate of A4 C from a global analysis of these observables except A was 



A was used in [ |2S|] , and an independent 
3|. In Fig. | we show the allowed core 



performed in [ |2§|] and re-analysed in [ 
analysis using all four observables in 
mass excess 5Ai c and envelope helium mass fraction Y env from the analyses [ 
Q |30|1 ; references to the original observations are found in these papers. 

Figure [5| suggests that, within the given uncertainties, the different observa- 
tions overlap at the standard core mass (5A4 C = 0) and at an envelope helium 
abundance of Y env which is compatible with the primordial helium abundance. 
Of course, the error bands do not have a simple interpretation because they com- 
bine observational and estimated systematic errors, which involve some subjective 
judgement by the authors. The difference between the two panels of Fig. || gives 
one a sense of how sensitive the conclusions are to these more arbitrary aspects 
of the analysis. As a nominal limit it appears safe to adopt |<5.M C | ^5 0.025 or 
\5M C \/ M. c ^ 5%; how much additional "safety-margin" one wishes to include is 
a somewhat arbitrary decision which is difficult to make objective in the sense of 
a statistical confidence level. 

In [ U it was shown that this limit can be translated into an approximate limit 
on the average anomalous energy-loss rate e x of a helium plasma, 

e x < 10 erg g" 1 s" 1 at T « 10 s K, p ss 2 x 10 5 g cm" 3 . (2) 

The density represents the approximate average of a red-giant core before helium 
ignition; the value at its center is about 10 6 g cm -3 . The main standard-model 




Figure 5: Allowed values for a core-mass excess at helium ignition 5M C and the 
envelope helium mass fraction Y env of evolved globular-cluster stars. Left after [ 
Q, right after [ 3C]. The observables are the brightness difference AM^ between 
the HB and the RGB tip, the RR Lyrae mass-to- light ratio A, their absolute 
brightness Mrr, and the number ratio R between HB and RGB stars. 
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neutrino emission process is plasmon decay 7 — > vv with a core average of about 
4 erg g" 1 s" 1 . Therefore, Eq. (Q) means that a new energy-loss channel must be 
less effective than a few times the standard neutrino losses. 



3.1.3 Helium-Burning Lifetime of Horizontal-Branch Stars 

We now turn to an energy-loss mechanism which becomes effective in a nonde- 
generate medium, i.e. we imagine that the core expansion after helium ignition 
"switches on" an energy-loss channel that was negligible on the RGB. Therefore, 
the pre-HB evolution is taken to be standard. As in the case of the Sun (Sec. pO| ) 
there will be little change in the HB stars' brightness, rather they will consume 
their nuclear fuel faster and thus begin to ascend the AGB sooner. The net 
observable effect is a reduction of the number of HB relative to RGB stars. 

From the measured HB/RGB number ratios in 15 globular clusters [31] and 
with plausible assumptions about the uncertainties of other parameters one con- 
cludes that the duration of helium burning agrees with stellar-evolution theory to 
within about 10% [ This implies that the new energy loss of the helium core 
should not exceed about 10% of its standard energy production rate. Therefore, 
the new energy-loss rate at average core conditions is constrained by [ 13 



, < 10 



erg g 



at T ss 0.7 x 10 8 K, 



P 



0.6 x 10 4 g cm" 3 . 



(3) 



This limit is slightly more restrictive than the often-quoted "red-giant bound," 
corresponding to e x < 100 erg g" 1 s" 1 at T = 10 8 K and p 



10 4 g cm 3 . 



It 



was based on the helium-burning lifetime of the "clump giants" in open clusters [ 
32]. They have fewer stars, leading to statistically less significant limits. The 
"clump giants" are the physical equivalent of HB stars, except that they occupy 
a common location at the base of the RGB, the "red-giant clump." 



3.1.4 Applications 

After the energy-loss argument has been condensed into the simple criteria of 
Eqs. (||) and (||) it can be applied almost mechanically to a variety of cases. 
The main task is to identify the dominant emission process for the new particles 
and to calculate the energy-loss rate e x for a helium plasma at the conditions 
specified in Eqs. @ or @. The most important limits will be discussed in the 
context of specific particle-physics hypotheses in Sees. |5|-f7|. Here we just mention 
that these and similar arguments were used to constrain neutrino electromagnetic 

|, ||, ||, ||, |7| , axions [ ||, ||, ||, g H g|, H> II rH 
] , paraphotons [ ^3l , the photo production cross section 



properties 
HP 0,14 




on 



|2 

He of new bosons 



54, p5|, the Yukawa couplings of new bosons to baryons 



or electrons [ 56, 57], and supersymmetric particles [ |5g, 59, 30] 

One may also calculate numerical evolution sequences including new energy 
18| , 2£ , |37|, j2], OT]. Comparing the results from such studies with what 
and (|3|) reveals that, in view of the overall theoretical and 



losses 

one finds from Eqs 

observational uncertainties, it is indeed enough to use these simple criteria 



3.2 White Dwarfs 

White dwarfs are another case where astronomical observations provide useful 
limits on new stellar energy losses. These compact objects are the remnants of 
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Figure 6: Observed white-dwarf luminosity function [ 33, 64]. The dotted line 



represents Mestel's cooling law with a constant birthrate of 10 -3 pc~ 3 Gyr -1 . 
The dashed line is from the cooling curve of a 0.6 M.q white dwarf which includes 
neutrino losses ] |65|], assuming the same constant birthrate. 



stars with initial masses of up to several M@ [ 27 , |6^ ]. For low-mass progen- 



itors the evolution proceeds as described in Sec. |3.1.1| . When they ascend the 
asymptotic giant branch they eventually shed most of their envelope mass. The 
degenerate carbon-oxygen core, having reached something like 0.6 M®, never 
ignites. Its subsequent evolution is simply one of cooling, first dominated by 
neutrino losses throughout its volume, later by surface photon emission. 

The cooling speed can be observationally infered from the "luminosity func- 
tion," i.e. the white-dwarf number density per brightness interval. As white 
dwarfs are intrinsically dim they are observed only in the solar neighborhood, 
out to perhaps 100 pc (1 pc = 3.26 lyr) which is far less than the thickness of the 
galactic disk. The measured luminosity function (Fig. |6|) reveals that there are 
few bright white dwarfs and many faint ones. The dotted line represents Mes- 



tel's cooling law [ 62, pq], an analytic treatment based on surface photon cooling. 



The observed luminosity function dips at the bright end, a behavior ascribed to 
neutrino emission which quickly "switches off" as the star cools. 

The luminosity function drops sharply at the faint end. Even the oldest white 
dwarfs have not yet cooled any further, implying that they were born 8-12 Gyr 
ago, in good agreement with the estimated age of the galaxy. Therefore, a novel 
cooling agent cannot be much more effective than the surface photon emission. 
This conclusion also follows from the agreement between the implied birthrate 
with independent estimates. The shape of the luminosity function can be de- 
formed for an appropriate temperature dependence of the particle emission rate, 
e.g. enhancing the "neutrino dip" at the bright end. Finally, white dwarfs in 
a certain range of surface temperatures are pulsationally unstable and are then 
called ZZ Ceti stars. The pulsation period of a few minutes depends on the lu- 
minosity, the period decrease thus on the cooling speed. For G117-B15A the 



period change was measured [ 67, nq], implying a somewhat large cooling rate. 



While this discrepancy may be worrisome, probably these measurements should 
be taken as an approximate confirmation of the predicted white-dwarf cooling 
speed. 
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White dwarfs were used to constrain the axion-electron coupling [ 69, 70, 71 
72, 73, 74]. It was also noted that the somewhat large period decrease of Gil 7- 
B15A could be ascribed to axion cooling [ |7q]. Finally, a limit on the neutrino 
magnetic dipole moment was derived [ [fq| . A detailed review of these limits 
is provided in [ Q; they are somewhat weaker than those from globular-cluster 
stars, but on the same general level. Therefore, white-dwarf cooling essentially 
corroborates some of the globular-cluster limits, but does not improve on them. 

3.3 Old Neutron Stars 

Neutron stars are the compact remnants of stars with initial masses beyond about 
8A4q. After their formation in a core-collapse supernova (Sec. ||) they evolve 
by cooling, a process that speeds up by a new energy-loss channel. Neutron-star 
cooling can now be observed by satellite-borne x-ray measurements of the thermal 



surface emission of several old pulsars — a recent review is [ 76]. 



Limits on axions were derived in [ 77, 78, 79], on neutrino magnetic dipole 



moments in f SO] . These bounds are much weaker than those from SN 1987A 



or globular clusters. Turning this around, anomalous cooling effects by particle 
emission is probably not important in old neutron stars, leaving them as labora- 
tories for many of the other uncertain bits of input physics such as the existence 
of new phases of nuclear matter [|8|, ||], |76|, [H], 32 1. If a neutron star converts into 



a strange-matter star an axion burst emerges [83], but for now this effect has 



not provided new empirical information on axion properties. 

4 SUPERNOVAE 

4-1 SN 1 981 A Neutrino Observations 

When the explosion of the star Sanduleak —69 202 was detected on 23 Febru- 
ary 1987 in the Large Magellanic Cloud, a satellite galaxy of our Milky Way at a 
distance of about 50 kpc (165,000 lyr), it became possible for the first time to mea- 
sure the neutrino emission from a nascent neutron star, turning this supernova 
(SN 1987A) into one of the most important stellar particle-physics laboratories [ 

A type II supernova explosion [ 87, |88], ^9], [x], 91, 92] is physically the implosion 



of an evolved massive star (A4 ^ 8 Mq) which has become an "onion-skin struc- 
ture" with several burning shells surrounding a degenerate iron core. It cannot 
gain further energy by fusion so that it becomes unstable when it has reached the 
limiting mass (Chandrasekhar mass) of 1-2 Mq that can be supported by elec- 
tron degeneracy pressure. The ensuing collapse is intercepted when the equation 
of state stiffens at around nuclear density (3 x 10 14 g cm -3 ), corresponding to a 
core size of a few tens of kilometers. At temperatures of tens of MeV this compact 
object is opaque to neutrinos. The gravitational binding energy of the newborn 
neutron star ("proto neutron star") of about 3 x 10 53 erg is thus radiated over 
several seconds from the "neutrino sphere." Crudely put, the collapsed SN core 
cools by thermal neutrino emission from its surface. 

The neutrino signal from SN 1987A (Fig. |?]) was observed by the v e p — > ne + 
reaction in several detectors [ ^6[. The number of events, their energies, and the 
distribution over several seconds corresponds well to theoretical expectations and 
thus has been taken as a confirmation of the standard picture that a compact 
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Figure 7: SN 1987A neutrino observations at Kamiokande [ p3|| , 1MB [ and 
Baksan [ |9f|]. The energies refer to the secondary positrons from the reaction 
v e p — > ne + . In the shaded area the trigger efficiency is less than 30%. The clocks 
have unknown relative offsets; in each case the first event was shifted to t = 0. 
In Kamiokande, the event marked as an open circle is attributed to background. 



remnant formed which emitted its energy by quasi-thermal neutrino emission. 



Detailed statistical analyses of the data were performed in [ 96, p7| . 

The signal does show a number of "anomalies." The average v e energies infered 
from the Irvine-Michigan-Brookhaven (1MB) and Kamiokande observations are 
quite different [ 9£j. The large time gap of 7.3 s between the first 8 and the last 
3 Kamiokande events looks worrisome [ |100| ] . The distribution of the final-state 
positrons from the i> e p — > ne + capture reaction should be isotropic, but is found 



to be significantly peaked away from the direction of the SN [ 55, 101, 102]. In 
the absence of other explanations, these features have been blamed on statistical 
fluctuations in the sparse data. 



4-2 Signal Dispersion 

A dispersion of the neutrino burst can be caused by a time-of-flight delay from a 
nonvanishing neutrino mass [ 103|. The arrival time from SN 1987A at a distance 
D would be delayed by 



At = 2.57 s 



D 



50 kpc 



10 MeV\ 2 



EL 



m,, 



10 eV 



(4) 



As the v e were registered within a few seconds and had energies in the 10 MeV 
range, m Ve is limited to less than around 10 eV. Detailed analyses reveal that 
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the pulse duration is consistently explained by the intrinsic SN cooling time and 
that m„„ < 20 eV is implied as something like a 95% CL limit [ 96], |104 |. 



The apparent absence of a time-of-flight dispersion effect of the v e burst was 
also used to constrain a "millicharge" of these particles (they would be deflected 
in the galactic magnetic field) |105[| , a quantum field theory with a fundamental 
length scale [ |106| ], and deviations from the Lorentzian rule of adding velocities [ 
107]. Limits on new long-range forces acting on the neutrinos [ 108, 109, 111 



111 



112] seem to be invalidated in the most interesting case of a long-range leptonic 



113 



force by screening from the cosmic background neutrinos 

The SN 1987A observations confirm that the visual SN explosion occurs several 
hours after the core-collapse and thus after the neutrino burst. Again, there is 
no apparent time-of-flight delay of the relative arrival times between the neutrino 
burst and the onset of the optical light curve, allowing one to confirm the equality 
of the relativistic limiting velocity for these particle types to within 2 x 10 -9 [ 
|114| , |115| ]. Moreover, the Shapiro time delay in the gravitational field of the 
galaxy of neutrinos agrees with that of photons to within about 4 x 10 -3 
constraining certain alternative theories of gravity [ 117, 118 ], 



4-3 Energy-Loss Argument 

The late events in Kamiokande and 1MB reveal that the signal duration was 
not anomalously short. Very weakly interacting particles would freely stream 
from the inner core, removing energy which otherwise would power the late- 
time neutrino signal. Therefore, its observed duration can be taken as evidence 
against such novel cooling effects. This argument has been advanced to constrain 
axion-nucleon couplings [ [HI [l20|, [l2l] , |l22] , p|, |24|, |||, [l26|, |i~27l majorons [ 
123, [1291 |l30l EML ^,1m,T3J, supersymmetric particles Hp5l, PI [jM 




higher dimensions 

neutrinos interacting by a Dirac mass term 
|152j , 153, 154, 155], mixed with active neutrinos 
right-handed currents [ 12C, 158, 15£, |l6 
HIT 



and graviton emission in quantum-gravity theories with 
It has also been used to constrain right-handed 

pc|, pBl pgl E37L pi, p|, pel, [151 



163, 164, 165, 166], or an electric form factor [ 167, 168 



156 , 157], interacting through 
161], a magnetic dipole moment [ 
Many of these 



results will be reviewed in Sees. ||-[7] in the context of specific particle-physics 
hypotheses. 

Here we illustrate the general argument with axions (Sec. ||) which are produced 
by nucleon bremsstrahlung AW — > NNa so that the energy-loss rate depends on 
the axion-nucleon Yukawa coupling g a ^ . In Fig. || we show the expected neutrino- 
signal duration as a function of g a N- With increasing g a N, corresponding to an 
increasing energy-loss rate, the signal duration drops sharply. For a sufficiently 
large g a N, however, axions no longer escape freely; they are trapped and thermally 
emitted from the "axion sphere" at unit optical depth. Beyond some coupling 
strength axions are less important than neutrinos and cannot be excluded. 

However, particles which are on the "strong interaction" side of this argument 
need not be allowed. They could be important for the energy-transfer during 
the infall phase and they could produce events in the neutrino detectors. For 
example, "strongly coupled" axions in a large range of g a N are actually excluded 
because they would have produced too many events by their absorption on 16 [ 

Likewise, particles on the free-streaming side can cause excess events in the 
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neutrino detectors. For example, right-handed neutrinos escaping from the inner 
core could become "visible" by decaying into left-handed states [ |170| ] or by spin- 
precessing in the galactic magnetic field if they have a dipole moment. 

Returning to the general argument, one can estimate a limit on the energy- 
loss rate on the free-streaming side by the simple criterion that the new channel 
should be less effective than the standard neutrino losses, corresponding to [ P§] 



< 10 



19 



erg g 



at p = 3 x 10 14 g cm 



30 MeV. 



(5) 



The density is the core average, the temperature an average during the first few 
seconds. Some authors find higher temperatures, but for a conservative limit it 
is preferable to stick to a value at the lower end of the plausible range. At these 
conditions the nucleons are partially degenerate while the electrons are highly 
degenerate. Several detailed numerical studies reveal that this simple criterion 
corresponds to approximately halving the neutrino signal duration [ Q 



121 | 



A simple analytic treatment is far more difficult on the trapping side; see 
for an example in the context of axions. 

The SN 1987A energy-loss argument tends to be most powerful at constraining 
new particle interactions with nucleons. Therefore, it is necessary to calculate 
the interaction rate with a hot and dense nuclear medium that is dominated 
by many-body effects. Besides the sparse data, the theoretical treatment of the 
emission rate is the most problematic aspect of this entire method. 



4-4 Radiative Neutrino Decays 

If neutrinos have masses one expects that the heavier ones are unstable and decay 
radiatively as v — > 1/7. SN 1987A is thought to have emitted similar fluxes of 
neutrinos and antineutrinos of all flavors so that one would have expected a burst 
of 7-rays in coincidence with the neutrinos. No excess counts were observed in the 
gamma-ray spectrometer (GRS) on the solar maximum mission (SMM) satellite [ 



171 



173, 174, 175 



171, 172], leading to restrictive limits on neutrino decays 
The GRS happened to go into calibration mode about 223 s after the neutrino 
burst, but for low- mass neutrinos (m u ^ 40 eV) the entire 7-ray burst would 
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Figure 8: Relative duration of SN neutrino cooling as a function of the axion- 
nucleon coupling. Freely streaming axions are emitted from the entire core vol- 
ume, trapped ones from the "axion sphere." The solid line follows from the 



numerical calculations [ 124, 125|; the dotted line is an arbitrary continuation. 
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have been captured, leading to a radiative decay limit of [ || 

r 7 /m„ > 0.8 x 10 15 s/eV. (6) 

For higher-mass neutrinos the photon burst would have been stretched beyond 
the GRS window. The first 7-rays from decays near the SN would arrive in 
coincidence with the v e burst, but the 7-burst duration would be given by some- 
thing like Eq. (Q). As a further complication, such higher-mass neutrinos violate 
the cosmological mass limit unless they decay sufficiently fast and thus nonra- 
diatively. Put another way, one must simultaneously worry about radiative and 
nonradiative decay channels — a detailed discussion is in [ |j . 

Comparable limits in the higher-mass range were also derived from 7-ray data 
of the Pioneer Venus Orbiter (PVO) [ |176f| . For m u > 0.1 MeV, decay photons 
still arrive years after SN 1987A. In 1991 the COMPTEL instrument aboard the 
Compton Gamma Ray Observatory looked at the SN 1987A remnant for about 
0.68 x 10 6 s, providing the most restrictive limits in this mass range [ |177| , |178f| . 

For m v <; 2m e ~ 1.2 MeV, which is only possible for v T with an experimental 
mass limit of about 18 MeV, the dominant radiative decay channel is v T — » 
v e e + e~ . From SN 1987A one would still expect 7-rays from the bremsstrahlung 



process v T — > i , e e + e~j, leading to interesting limits [ 172, 176, 179, |180| , 181]. 

The decay positrons from past galactic SNe would be trapped by the galactic 
magnetic fields and thus linger for up to 10 5 yr. Independently of SN 1987A, 
measurements of the galactic positron flux thus provide limits on neutrino decays 
with final-state positrons [ [|, 182]. 



4-5 Explosion Energetics 

The standard scenario of a type II SN explosion has it that a shock wave forms 
near the edge of the core when its collapse halts at nuclear density and that this 
shock wave ejects the mantle of the progenitor star. However, in typical numerical 
calculations the shock wave stalls so that this "prompt explosion" scenario does 
not seem to work. In the "delayed explosion" picture the shock wave is revived 
by neutrino heating, perhaps in conjunction with convection, but even then it 
appears difficult to obtain a successful or sufficiently energetic explosion. 

Therefore, one may speculate that nonstandard modes of energy transfer play 
an important role. An example is Dirac neutrinos with a magnetic dipole moment 
of order 10~ 12 ub (Bohr magnetons). The right-handed (sterile) components 
would arise in the deep inner core by helicity-flipping collisions and escape. They 
precess back into interacting states in the large magnetic fields outside the SN 
c^ore and heat the shock region; their interaction cross section would be relatively 



large because of their large inner-core energies [ [L83| , |184| , |185| , [L86|, 

Certainly it is important not to deposit too much energy in the mantle and 
envelope of the star. 99% of the gravitational binding energy of the neutron star 
goes into neutrinos, about 1% into the kinetic energy of the explosion, and about 
0.01% into the optical supernova. Therefore, neutrinos or other particles emitted 
from the core must not decay radiatively within the progenitor's envelope radius 



of about 100 s or else too much energy lights up [ 189, 09 



18 



Georg G. Raffelt 



4-6 Neutrino Spectra and Neutrino Oscillations 



Neutrino oscillations can have several interesting ramifications in the context 
of SN physics because the temporal and spectral characteristics of the emission 
process depend on the neutrino flavor [ 90, 91, |9^, 203]. The simplest case is that 
of the "prompt u e burst" which represents the deleptonization of the outer core 
layers at about 100 ms after bounce when the shock wave breaks through the edge 
of the collapsed iron core. This "deleptonization burst" propagates through the 
mantle and envelope of the progenitor star so that resonant oscillations take place 
for a large range of mixing parameters between v e and some other flavor, notably 



for most of those values where the MSW effect operates in the Sun [ 191, 192, 193 
19l , IH, IH, |l97|, HH, |l9|, |(X|. In a water Cherenkov detector this burst is 



visible as v e -e scattering, which is forward peaked, but one would have expected 
only a fraction of an event from SN 1987A. The first event in Kamiokande may 
be attributed to this signal, but this interpretation is statistically insignificant. 

During the next few hundred milliseconds the shock wave stalls at a few hun- 
dred kilometers above the core and needs rejuvenating. The efficiency of neutrino 
heating can be increased by resonant flavor oscillations which swap the v e flux 
with, say, the v T one. Therefore, what passes through the shock wave as a v e 
was born as a u T at the proto neutron star surface. It has on average higher 
energies and thus is more effective at transfering energy. In Fig. [9| the shaded 
range of mixing parameters is where supernovae are helped to explode, assuming 
a "normal" neutrino mass spectrum with m Ve < m Ur [ 201]. Below the shaded 
region the resonant oscillations take place beyond the shock wave and thus do 
not affect the explosion. 
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Figure 9: Mass difference and mixing between v e and z/« or u T where a spectral 
swap would occur to help explode supernovae, schematically after f |201 
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215, 216 



217 



The logic of this scenario depends on deviations from strictly thermal neutrino 
emission at some blackbody "neutrino sphere." The neutrino cross sections are 
very energy dependent and different for different flavors so that the concept of 
a neutrino sphere is rather crude — the spectra are neither thermal nor equal 
for the different flavors [ 202, 202]. The dominant opacity source for v e is the 
process v e + n — > p + e~ , for v e it is v e + p — > n + e + , while for v^ T and z? M)T 
it is neutral-current scattering on nucleons. Therefore, unit optical depth is at 
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the largest radius (and lowest medium temperature) for v e , and deepest (highest 
temperature) for f „ T and T . In typical calculations one finds a hierarchy 
(E Ue ) : (E Pe ) : (£ othere ) « | : 1: | with (£ Pe ) = 14-17 MeV [||. The SN 1987A 
observations imply a somewhat lower range of (Ep e ) ~ 7-14 MeV [ 96, J99|. 

It should be noted that, pending a more detailed numerical confirmation [ 
204], the difference between the v e and v^ T or v^ T average energies appears to 
be smaller than commonly assumed [ 126, 205, 206, 207], but there is no doubt 



that the v e spectrum is softer than the others. Still, the quantitative import of 
flavor oscillations depends on details of the neutrino spectra formation process in 
those SN core layers where the diffusion approximation for the neutrino transport 
is no longer valid, yet neutrinos are still trapped. 

A few seconds after core bounce the shock wave has long since taken off, leaving 
behind a relatively dilute "hot bubble" above the neutron-star surface. This 
region is one suspected site for r-process heavy-element synthesis, which requires 



a neutron-rich environment [ 208 , 209 , 210 , 211 , 212 , |213| , 214 ]. The neutron- 
to-proton ratio, which is governed by the (5 reactions u e + n — > p + e~ and 
v e + p — ► n + e + , is shifted to a neutron-rich phase if (E Ue ) < (Ep e } as for 
standard neutrino spectra. Resonant oscillations can again swap the v e flux with 
another one, inverting this hierarchy of energies. In the hatched range of mixing 
parameters shown in Fig. || the r-process would be disturbed [ 215, [216| , 217 , 218 1 . 



On the other hand, v e — * v s oscillations into a sterile neutrino could actually help 
the r-process by removing some of the neutron-stealing v e [ 21S, |236[| . 



A large body of recent literature was devoted to explaining the large kick 
velocities of the observed radio pulsars as a "neutrino rocket effect." The required 
few-percent anisotropy of the SN neutrino emission was attributed to an intricate 
interplay between the magnetic-field induced neutrino dispersion relation and 

223] , [224], [225], |226|, p^. However, due to 



resonant oscillations \ 220, 221 



228 1 so that the pulsar 



a conceptual error the effect was vastly overestimated 
kicks do not seem to be related to neutrino oscillations in any obvious way. 

If the mixing angle between v e and some other flavor is large, the v e flux from 
a SN contains a significant fraction of oscillated states that were born as 
or v T and thus should have higher average energies. The measured SN 1987A 
event energies are already somewhat low, a problem so strongly exacerbated by 
oscillations that a large-mixing-angle solution of the solar neutrino deficit poses 
a problem 



39], 104, 229]. This conclusion, however, depends on the standard 
predictions for the average neutrino energies which may not hold up to closer 
scrutiny as mentioned above. 



5 LIMITS ON NEUTRINO PROPERTIES 



5.1 Masses and Mixing 

Astrophysics and cosmology play a fundamental role for neutrino physics as the 
properties of stars and the universe at large provide some of the most restrictive 
limits on nonstandard properties of these elusive particles. Therefore, it behoves 
us to summarize what the astrophysical arguments introduced in the previous 
sections teach us about neutrinos. 

Unfortunately, stars do not tell us very much about neutrino masses, the holy 



grail of neutrino physics. The current discourse [ 23C, 231, 232] centers on the 
interpretation of the solar [ u\ and atmospheric [ 233] neutrino anomalies and 
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the LSND experiment [ 234, 235] which all provide very suggestive evidence for 
neutrino oscillations. Solar neutrinos imply a Am„ of about 10 -5 eV 2 (MSW solu- 
tions) or 10" 10 eV 2 (vacuum oscillations), atmospheric neutrinos 10 3 ~10 2 eV , 
and the LSND experiment 0.3-8 eV 2 . Taken together, these results require a 
fourth flavor, a sterile neutrino, which is perhaps the most spectacular implica- 
tion of these experiments, but of course also the least secure. 

Core-collapse supernovae appear to be the only case in stellar astrophysics, 
apart from the solar neutrino flux, where neutrino oscillations can be important. 
However, Fig. [9| reveals that the experimentally favored mass differences negate a 
role of neutrino oscillations for the explosion mechanism or r-process nucleosyn- 
thesis, except perhaps when sterile neutrinos exist 
the interpretation of the SN 1987A signal 
galactic SN 



219, 23§|. Oscillations affect 



237 



However, as discussed in Sec. 4.6 



104] , HU and_that of a future 
the main chal- 



lenge is to develop a quantitatively more accurate understanding of supernovae 
as neutrino sources before relying on relatively fine points of the neutrino spectral 
characteristics to learn about neutrino mixing parameters. 

Oscillation experiments reveal only mass differences so that one still needs 
to worry about the absolute neutrino mass scale. The absence of anomalous 



SN 1987A signal dispersion (Sec. 4.2) gives a limit [ |96|, 104] 



£ 20 eV, 



(7) 



somewhat weaker than current laboratory bounds. A high-statistics observation 
of a galactic SN by a detector like Superkamiokande could improve this limit 
to about 3 eV by using the fast rise-time of the neutrino burst as a measure of 
dispersion effects [ 240]. If the neutrino mass differences are indeed very small, 
this limit carries over to the other flavors. One can derive an independent mass 
limit on v „ and v T in the range of a few 10 eV if one identifies a neutral-current 



signature in a water Cherenkov detector [ 241, 242, 243, 244, 245], or if a future 
neutral-current detector provides an additional measurement [ 24£ , |247| ] . 

The SN 1987A energy-loss argument (S ec. |4,3j ) provides a limit on a neutrino 
Dirac mass of [ |, |12(], [L45|, |M§, [L47|, [L4§ pjl 

m,, (Dirac) < 30 keV. (8) 

It is based on the idea that trapped Dirac neutrinos produce their sterile com- 
ponent with a probability of about (m u /2E u ) 2 in collisions and thus feed energy 
into an invisible channel. This result was important in the discourse on Simpson's 
17 keV neutrino which is now only of historical interest 



245 



5.2 Dipole and Transition Moments 
5.2.1 Electromagnetic Form Factors 

Neutrino electromagnetic interactions would provide for a great variety of astro- 
physical implications. In the vacuum, the most general neutrino interaction with 
the electromagnetic field is [ [249| , [250| 

Ant = -Fifo^A" - Gi^vys^ V* - \^^{F 2 + G 275 )V^, (9) 

where tp is the neutrino field, A 1 the electromagnetic vector potential, and F^ v 
the field-strength tensor. The form factors are functions of Q 2 with Q the energy- 
momentum transfer. In the Q 2 — > limit F\ is the electric charge, G\ an anapole 
moment, Fi a magnetic, and Gi an electric dipole moment. 
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If neutrinos are electrically strictly neutral, corresponding to -F\(0) = 0, they 
still have a charge radius, usually defined as (r 2 ) = 6dFi(Q 2 )/edQ 2 \Q2 =0 . This 
form factor provides for a contact interaction, not for a long-range force, and 
as such modifies processes with Z° exchange 



251, 



254 |255| - As as 



trophysics provides no precision test for the effective strength of neutral-current 
interactions, this form factor is best probed in laboratory experiments [ 256 1. 
Likewise, the anapole interaction vanishes in the Q 2 — > limit and thus repre- 
sents a modification to the standard neutral-current interaction, with no apparent 
astrophysical consequences. 

The most interesting possibility are magnetic and electric dipole and transition 
moments. If the standard model is extended to include neutrino Dirac masses, the 
magnetic dipole moment is fi u = 3.20 x 10~ 19 {j,b m u /eV where (ib = e/2m e is the 
Bohr magneton 
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50j. An electric dipole moment e„ violates CP, and both 
are forbidden for Majorana neutrinos. Including flavor mixing implies electric and 
magnetic transition moments for both Dirac and Majorana neutrinos, but they 
are even smaller due to GIM cancellation. These values are far too small to be of 
any experimental or astrophysical interest. Significant neutrino electromagnetic 
form factors require a more radical extension of the standard model, for example 
the existence of right-handed currents. 



5.2.2 Plasmon Decay in Stars 

Dipole or transition moments allow for several interesting processes (Fig. |l0|). For 
the purpose of deriving limits, the most important case is 7 — > vv which is kine- 
matically possible in a plasma because the photon acquires a dispersion relation 
which roughly amounts to an effective mass. Even without anomalous couplings, 
the plasmon decay proceeds because the charged particles of the medium induce 
an effective neutrino-photon interaction. Put another way, even standard neu- 



trinos have nonvanishing electromagnetic form factors in a medium [ 257, 258 1. 
The standard plasma process [ 259, [260| , 261 1 dominates the neutrino production 



in white dwarfs or the cores of globular-cluster red giants. 
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Figure 10: Processes with neutrino electromagnetic dipole or transition moments. 
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The plasma process was first used in [ |262| ] to constrain neutrino electromag- 
netic couplings. Numerical implementations of the nonstandard rates in stellar- 
evolution calculations are [ ^9], 35, [37], The helium- ignition argument in 
globular clusters (Sec. 3.1. 2D , equivalent to Eq. (||), implies a limit [[|, |2^, ^g, |37| 



iiv % 3 x 10 



-12 



MB, 



(10) 



applicable to magnetic and electric dipole and transition moments for Dirac and 
Majorana neutrinos. Of course, the final-state neutrinos must be lighter than the 
photon plasma mass which is around 10 keV for the relevant conditions. 

The corresponding laboratory limits are much weaker [ ||. The most restrictive 
bound is \i Ve < 1.8 x 10~ 10 /xb at 90% CL from a measurement of the i> e -e- 
scattering cross section involving reactor sources. A significant improvement 
should become possible with the MUNU experiment [ 265], but it is unlikely that 



the globular-cluster limit can be reached anytime soon. 



5.2.3 Radiative Decay 

A neutrino mass eigenstate v\ may decay to another one Vj by the emission of a 
photon, where the only contributing form factors are the magnetic and electric 



transition moments. The inverse radiative lifetime is found to be [ 249, 250 1 



m A 
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where /i« and e« are the transition moments while |/x e ff| 
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my 
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. Radia- 



tive neutrino decays have been constrained from the absence of decay photons of 



reactor i> e fluxes [ 264 1, the solar v e flux [ [265| , 266 1, an d the SN 1987A neutrino 
burst [ 171, 172, 173, 174, 175 ]. For m v = m; ^> rrij these limits can be expressed 
as 
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Sun (z; e ), 

SN 1987A (all flavors), 
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(12) 



In this form the SN 1987A limit applies for m v % 40 eV as explained in Sec. 4.4 . 
The decay of cosmic background neutrinos would contribute to the diffuse photon 
backgrounds, excluding the shaded areas in Fig. 11 
delineated by the dashed line 



They are approximately 
corresponding to the bottom line in Eq. (^%). 



More restrictive limits obtain for certain masses above 3 eV from the absence of 
emission features from several galaxy clusters 



270, 271, B8 



For low-mass neutrinos the phase-space factor in Eq. (|ll]) is so punishing 
that the globular-cluster limit is the most restrictive one for m v below a few 
eV. This is precisely the mass range which today appears favored from neutrino 
oscillation experiments. Turning this around, the globular-cluster limit implies 
that radiative decays of low-mass neutrinos do not seem to have observable con- 
sequences. 

For masses above about 30 eV one must invoke fast invisible decays in order 
to avoid a conflict with the cosmological mass limit. In this case radiative decay 
limits involve the total lifetime as another parameter; the SN 1987A limits have 
been interpreted in this sense in [0, 174, 176, |178||. 
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Figure 11: Astro-physical limits on neutrino dipole moments. The light-shaded 



background-radiation limits are from [ 267 1, the dark-shaded ones from [ 265 , 269 1, 
the dashed line is the approximation formula in Eq. (12), bottom line. 



5.2.4 Cherenkov Effect 

Another form of "radiative decay" is the Cherenkov effect v — ► v + 7, which 
involves the same initial- and final-state neutrino. This process is kinematically 
allowed for photons with uj 2 — k 2 < 0, which obtains in certain media or in 
external magnetic fields. The neutrino may have an anomalous dipole moment, 
but there is also a standard-model photon coupling induced by the medium or 
the external field. Thus far it does not look as if the neutrino Cherenkov effect 



has any strong astrophysical significance (see [ 272 ] for a review of the literature) . 
5.2.5 Spin-Flip Scattering 

The magnetic or electric dipole interaction couples neutrino fields of opposite chi- 
rality. In the relativistic limit this implies that a neutrino flips its helicity in an 
"electromagnetic collision," which in the Dirac case produces the sterile compo- 
nent. The active states are trapped in a SN core so that spin-flip collisions open 
an energy-loss channel in the form of sterile states. Conversely, the SN 1987A 
energy-loss argument (Sec. |4.3|) allows one to derive a limit [ |163|, |l6f 



fi u (Dirac) < 3 x 10 12 /xb, 



(13) 



for both electric and magnetic dipole and transition moments. It is the same as 
the globular-cluster limit Eq. (|l0|), which however includes the Major ana case. 

Spin-flip collisions would also populate the sterile Dirac components in the early 
universe and thus increase the effective number of thermally excited neutrino 
degrees of freedom. Full thermal equilibrium is attained for fi v (Dirac) ^ 60 x 
10~ 12 /jb [ H HH- In view of the SN 1987A and globular-cluster limits this 



result assures us that big-bang nucleosynthesis remains undisturbed. 
5.2.6 Spin and Spin-Flavor Precession 

Neutrinos with magnetic or electric dipole moments spin-precess in external mag- 



netic fields [ 274, 275]. For example, solar neutrinos can precess into sterile and 
thus undetectable states in the Sun's magnetic field [ 276, 277, 278]. The same 
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for SN neutrinos in the galactic magnetic field where an important effect obtains 
for [i v > 10~ 12 ^b- Moreover, the high-energy sterile states emitted by spin-flip 
collisions from the inner SN core could precess back into active ones and cause 
events with anomalously high energies in SN neutrino detectors, an effect which 
probably requires /^(Dirac) < lCP 12 /ie from the SN 1987A signal [ [163], |279| . 
For the same general [i v magnitude one may expect an anomalous rate of energy 
transfer to the shock wave in a SN, helping with the explosion (Sec. 



In a medium the refractive energy shift for active neutrinos relative to sterile 
ones creates a barrier to the spin precession [ 281, 282]. The mass difference has 
the same effect if the precession is between different flavors through a transition 
moment [ 28C ] . However, the mass and refractive terms may cancel, leading to res- 
onant spin- flavor oscillations in the spirit of the MSW effect [ 283 , 284 , 285 , 286 1. 



This mechanism can explain all solar neutrino data [ 287 , 288 1 , but requires rather 
large toroidal magnetic fields in the Sun since the neutrino magnetic (transition) 
moments have to obey the globular-cluster limit of Eq. (|To|). For Majorana neu- 
trinos, the spin- flavor precession amounts to transitions between neutrinos and 
antineutrinos so that the observation of antineutrinos from the Sun would be a 
diagnostic for this effect 



28g , |29C| , M - 



Large magnetic fields exist in SN cores so that spin-flavor precession could play 
an important role there, with possible consequences for the explosion mechanism, 
r-process nucleosynthesis, or the measurable neutrino signal [ |292| , 293, 294, 295, 
|296j ] . The downside of this richness of phenomena is that there are so many un- 
known parameters (electromagnetic neutrino properties, masses, mixing angles) 
as well as the unknown magnetic field strength and distribution that it is difficult 
to come up with reliable limits or requirements on neutrino properties. The SN 
phenomenon is probably too complicated to serve as a laboratory to pin down 
electromagnetic neutrino properties, but it clearly is an environment where these 
properties could have far-reaching consequences. 



5.3 Millicharged Particles 



It is conceivable that neutrinos carry small electric charges if charge conservation 

1 ML Eol. More- 



is not exact [ 297, 298] or if the families are not sequential 



over, new particles with small electric charges are motivated in certain models 
with a "mirror sector" and a slightly broken mirror symmetry [ 302]. Therefore, it 
is interesting to study the experimental, astrophysical, and cosmological bounds 

304]|355|, 



on "millicharged" particles [ t 

A model-independent v e charge limit arises from the absence of dispersion of 



307 ], 



the SN 1987A neutrino signal in the galactic magnetic field [ [|, |105| ] 



< 3 x l(T 17 e. 



(14) 



If charge conservation holds in neutron decay, e Ue ^ 
a limit for the neutron charge ofe n = (— 0.4 ± 1.1 
neutrality of matter which was found to be e p + e e 



3 x 10 21 e results, based on 
x 10" 21 e [ |H and on the 
= (0.8 ±0.8) x 10~ 21 e ]|309|. 



The measured Vn-e cross section implies e h 



< 



10- 



Generic millicharged particles (charge e x , mass m x ) could appear as virtual 
states and would thus modify the Lamb shift unless e x < O.llem^/MeV [ |304| . 
A number of limits follow from a host of previous accelerator experiments [ |304| 
and a recent dedicated search at SLAC [ 310] — see Fig. 12. 
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Millicharged particles are produced by the plasmon decay process and thus 
drain energy from stars. In globular clusters, the emission rate is almost the 
same for HB stars and red giants before helium ignition, in contrast with the 
magnetic-dipole case. Therefore, Eqs. (||) and (^) give an almost identical limit [ 



261] 



e x <2xW- lA e, (15) 

applicable for m x below a few keV. The SN 1987A energy-loss argument extends 
the exclusion range to about 10 MeV for 1(T 9 e < e x < 1(T 7 e [ f§ |i~67| . 

The usual big-bang nucleosynthesis (BBN) limit on the effective number of 
neutrino species N e Q provides another constraint. A millicharged neutrino is of 
Dirac nature so that its right-handed component adds one effective species. If the 
millicharged particles are not neutrinos, then depending on their spin N e Q may 
increase even more. If BBN excludes one extra species one finds e x ■ < 3 x 10 -9 e 
for m x <■ 1 MeV [ |304|| . More stringent limits apply in certain models where the 



millicharged particles are associated with a shadow sector [ |30T | . 



Further regions in Fig. 12 are excluded to avoid "overdosing" the universe by 
the new particles [ 304, 307 1. However, because their relic density depends on 
their annihilation cross section, it is necessary to specify a model. It is hard to 
imagine new particles which interact solely through their small electric charge! 



5.4 Nonstandard Weak Interactions 
5.4.1 Right-Handed Currents 

Right-handed (r.h.) weak interactions may exist on some level, e.g. in left-right 
symmetric models where the r.h. gauge bosons differ from the standard ones by 
their mass. In the low-energy limit relevant for stars one may account for the new 
couplings by a r.h. Fermi constant eG-p where e is a small dimensionless parameter. 
In left-right symmetric models one finds explicitly for charged-current processes 
e cc = C 2 + [ m (WL)/m(Wn)]' 2 where m(WL,R) are the l.h. and r.h. gauge boson 



masses and £ is the left-right mixing parameter [ 158]. 




Figure 12: Limits on the electric charge e x and mass m x of generic millicharged 
particles which may be neutrinos or new particles. In order to avoid overdosing 
the universe, additional model-dependent parameter regions are excluded. The 
big-bang nucleosynthesis (BBN) excluded region is larger in some models. 
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Assuming that neutrinos are Dirac particles, a SN core loses energy into r.h. 
states as an "invisible channel" by the process e + p —* n + v e ,R- The SN 1987A 
energy-loss argument (Sec. 4.3 ) then requires ecc ~ 10 -5 [ ||, 12C , |158] ]. Labora- 
tory experiments yield a weaker limit of order ecc ~ 3 x 10~ 2 [ 311], but do not 
depend on the assumed existence of r.h. neutrinos. 

For neutral currents the dominant emission process is NN — > NNvrvr which 
is subject to saturation effects as in the case of axion emission [ |126[] . One then 
finds £nc ~ 3 x 10~ 3 [p}|, somewhat less restrictive than the original limits of [ 
|12C|, tL58i; see also [ |159| , |160| , 161]. This bound is also somewhat less restrictive 



158| ; see also 

than £nc 10~ 3 found from big-bang nucleosynthesis 



5.4.2 Secret Neutrino Interactions and Majorons 

The neutrino-neutrino cross section is not known experimentally. It could be 
anomalously large if neutrino Majorana masses were to arise from a suitable 
majoron model [ 13131, 314, 315, B16fl. "Secret" neutrino-neutrino interactions 



were constrained by the fact that the SN 1987A neutrino signal was not depleted 
by collisions with cosmic background neutrinos [ 317j|. Supernova physics with 
majorons and SN 1987A limits were discussed in 

i~3l, |T|, Em 



320, 321 



li~2l , m H], |m|, m EH 



22] , |323f |. There is little doubt that majoron models will 
have an important impact on SN physics for neutrino- majoron Yukawa couplings 
in the 10 _6 -10~ 3 range. The existing literature, however, is too confusing for this 
author to come up with a clear synthesis of what SN physics implies for majoron 
models. 



5.4.3 Flavor-Changing Neutral Currents 

In certain models the neutrino neutral current has an effective flavor-changing 
component. Neutrinos propagating in matter then have medium-induced mixings 



and thus can oscillate even if they are strictly massless [ 324 , |325| . Naturally, this 
phenomenon can be important for the oscillation of solar 
and supernova 



330 , 331 1 neutrinos. 



6 AXIONS AND OTHER PSEUDOSCALARS 



6.1 Interaction Structure 

New spontaneously broken global symmetries imply the existence of Nambu- 
Goldstone bosons that are massless and as such present the most natural case (be- 
sides neutrinos) for using stars as particle-physics laboratories. Massless scalars 
would lead to new long-range forces (Sect. 0) so that we may focus here on pseu- 
doscalars. The most prominent example are axions which were proposed more 
than twenty years ago as a solution to the strong CP problem [ |332| , [333| , 334, 335 1; 
for reviews see [ 13361 337 ] and for the latest developments the proceedings of a 
topical conference [ 338]. We use axions as a generic example — it will be obvious 



how to extend the following results and discussions to other cases. 

Actually, axions are only "pseudo Nambu-Goldstone bosons" in that the spon- 
taneously broken chiral Peccei-Quinn symmetry Upq(1) is also explicitly broken, 
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providing these particles with a small mass 

m a = 0.60 eV — GeV . (16) 

J a 

Here, f a is the Peccei-Quinn scale, an energy scale which is related to the vacuum 
expectation value of the field that breaks C/pq(1). The properties of Nambu- 
Goldstone bosons are always related to such a scale which is the main quantity 
to be constrained by astrophysical arguments, while Eq. ( |T6| ) is specific to axions 
and allows one to express limits on f a in terms of m a . 

In order to calculate the axionic energy-loss rate from stellar plasmas one needs 
to specify the interaction with the medium constituents. The interaction with a 
fermion j (mass rrij) is generically 

Ant = wf- Vj-fysVjdna or - i Vj^ja, (17) 

where \Pj is the fermion and a the axion field and C~ is a model-dependent co- 
efficient of order unity. The combination g a j = Cjm,j/f a plays the role of a 
Yukawa coupling and a a j = g^/Air acts as an "axionic fine structure constant." 
The derivative form of the interaction is more fundamental in that it is invari- 
ant under a — > a + ag and thus respects the Nambu-Goldstone nature of these 
particles. The pseudoscalar form is usually equivalent, but one has to be care- 
ful when calculating processes where two Nambu-Goldstone bosons are attached 
to one fermion line, for example an axion and a pion attached to a nucleon [ 



12q , |339] |340| , |34|, |342 | 



The dimensionless couplings Cj depend on the detailed implementation of the 
Peccei-Quinn mechanism. Limiting our discussion to "invisible axion models" 
where f a is much larger than the scale of electroweak symmetry breaking, it is 
conventional to distinguish between models of the DFSZ type (Dine, Fischler, 
Srednicki [ |4|, Zhitnitskh [ |4|) and of the KSVZ type (Kim [ |345| , Shifman, 



Vainshtein, Zakharov [ [346f ) . In KSVZ models, axions have no tree-level couplings 
to the standard quarks or leptons, yet axions couple to nucleons by their generic 
mixing with the neutral pion. The latest analysis gives numerically [ 127| 



C p = -0.34, C n = 0.01 (18) 

with a statistical uncertainty of about ±0.04 and an estimated systematic un- 
certainty of roughly the same magnitude. The tree-level couplings to standard 
quarks and leptons in the DFSZ model depend on an angle (3 which measures the 
ratio of vacuum expectation values of two Higgs fields. One finds [ 127] 



C e = \ cos 2 p, C p = -0.07 - 0.46 cos 2 p, C n = -0.15 + 0.38 cos 2 p, (19) 

with similar uncertainties as in the KSVZ case. 

The CP-conserving interaction between photons and pseudoscalars is com- 
monly expressed in terms of an inverse energy scale g ai according to 

Ant = \g ai F^F^a = -g ai B ■ B a, (20) 

where F is the electromagnetic field-strength tensor and F its dual. For axions 

a E 

g a -y = ^f a c ~ti c ^ = n~ L92 ± a08 ' (21) 
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where E/N is the ratio of the electromagnetic and over color anomalies, a model- 
dependent ratio of small integers. In the DFSZ model or grand unified models 
one has E/N = 8/3, for which C 7 ~ 0.75, but one can also construct models 
with E/N = 2, which significantly reduces the axion-photon coupling [ 347]. The 



value of C 7 in a great variety of cases was reviewed in [ 348, [349| 



6.2 Limits on the Interaction Strength 
6.2.1 Photons 

The axion interaction with fermions or photons allows for numerous reactions 
which can produce axions in stars, which may imply limits on the axion cou- 
pling strength. Beginning with photons, pseudoscalars interact according to the 
Lagrangian of Eq. ( |20|) which allows for the decay a — > 27. In stellar plas- 
mas the photon-axion interaction also makes possible the Primakoff conversion 
7 <-> a in the electric fields of electrons and nuclei [ [3^] — see Fig. [l|. For low-mass 
pseudoscalars the emission rate was calculated for various degrees of electron 
74, 350 1, superseding an earlier calculation where screening 
@3|. 



degeneracy in [ 
effects had been ignored 

The helioseismological constraint on solar energy losses then leads to Eq. ([!]) 
as a bound on g a ^. Figure 12 shows this constraint ("Sun") in the context of 
other bounds; similar plots are found in [ |, § [351], |||, H|. For axions the 
relationship between g a -y and m a is indicated by the heavy solid line, assuming 
E/N = 8/3. 

One may also search directly for solar axions. One method ("helioscope") is 
to direct a dipole magnet toward the Sun, allowing solar axions to mutate into 




Figure 13: Limits to the axion-photon coupling g a ~ f as defined in Eq. (|20j). They 
apply to any pseudoscalar except for the "haloscope" search which assumes that 
these particles are the galactic dark matter; the dotted region marks the projected 
sensitivity range of the ongoing dark-matter axion searches. For higher masses 



than shown here the pertinent limits are reviewed in [ 352 1 . 
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x-rays by the inverse Primakoff process [ 355, 



J. A pilot experiment was not 
sensitive enough [ 357], but the exposure time was significantly increased in a new 
experiment in Tokyo where a dipole magnet was gimballed like a telescope so that 
it could follow the Sun [ |358 ], |359|1 . The resulting limit g ai ^6x 1CP 10 GeV -1 is 
more restrictive than Eq. Another helioscope project was begun in Novosi- 
birsk several years ago 



361 ], but its current status has not been reported for 
some time. An intruiging project (SATAN) at CERN would use a decommis- 
sioned LHC test magnet that could be mounted on a turning platform to achieve 
reasonable periods of alignment with the Sun ] |362| , This setup could begin to 
compete with the globular-cluster limit of Eq. (|22j). 

The axion-photon transition in a macroscopic magnetic field is analogous to 
neutrino oscillations and thus depends on the particle masses [ |360|] . For a large 
mass difference the transition is suppressed by the momentum mismatch of parti- 
cles with equal energies. Therefore, the Tokyo limit applies only for m a < 0.03 eV. 
In a next step one will fill the helioscope with a pressurized gas, giving the photon 
a dispersive mass to overcome the momentum mismatch. 

An alternative method is "Bragg diffraction," which uses uses the strong elec- 
tric field of a crystal lattice which has large Fourier components for the required 
momentum transfer 



363 



364] , |365[| . The experiment has been performed using 
Ge detectors which were originally built to search for neutrinoless double-beta 
decay and for WIMP dark matter; the crystal serves simultaneously as a Pri- 
makoff "transition agent" and as an x-ray detector. A first limit of the SOL AX 
Experiment [ |366[| of g a ~ f < 27 x 10~ 10 GeV -1 is not yet compatible with Eq. @ 
and thus not self-consistent. In the future one may reach this limit, but prospects 
to go much further appear dim [ 367]. 

The Primakoff conversion of stellar axions can also proceed in the magnetic 
fields of Sun spots or in the galactic magnetic field so that one might expect 
anomalous x- or 7-ray fluxes from the Sun [ |368j ] , the red supergiant Betelgeuse [ 
|369| 1, or SN 1987A [ g?(| Observations of SN 1987A yield g ai < 

10 -10 GeV -1 for nearly massless pseudoscalars with m a 



0.1 x 
A similar 



< 10~ 9 eV 

limit obtains from the isotropy of the cosmic x-ray background which would be 
modified by the conversion to axions in the galactic magnetic field [ |372| . Axion- 
photon conversion in the magnetic fields of stars, the galaxy, or the early universe 
were also studied in [ Hj| |7|, |7|, [377], |78|, |80| , but no additional 

limits emerged. 

The existence of massless pseudoscalars would cause a photon birefringence 
effect in pulsar magnetospheres, leading to a differential time delay between pho- 
tons of opposite helicity and thus to g ay < 0.5 x 10~ 10 GeV -1 [ |381[| . 

A laser beam in a laboratory magnetic field would also be subject to vacuum 
birefringence [ 382], adding to the QED Cotton-Mouton effect. First pilot ex- 
periments [ 353, |383f| did not reach the QED level. Two vastly improved current 
projects are expected to get there [ 384, |385| ] , but they will stay far away from the 
"axion line" in Fig. With a laser beam in a strong magnet one can also search 
for Primakoff axion production and subsequent back-conversion, but a pilot ex- 
periment naturally did not have the requisite sensitivity [ 38£]. The exclusion 
range of current laser experiments is schematically indicated in Fig. [T^. 

The most important limit on the photon coupling of pseudoscalars derives from 
the helium-burning lifetime of HB stars in globular clusters, i.e. from Eq. (0), 



9a-y 



< 0.6 x 10" 10 GeV 



(22) 
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For m a > 10 keV this limit quickly degrades as the emission is suppressed when 
the particle mass exceeds the stellar temperature. For a fixed temperature, the 
Primakoff energy-loss rate decreases with increasing density so that Eq. (|2|) im- 
plies a less restrictive constraint. Equation (^) was first stated in [Q, supersed- 
ing the slightly less restrictive but often-quoted "red-giant bound" of [ |l8| — see 
the discussion after Eq. (0). The axion relation Eq. (|2l|) leads to 



m a C<y < 0.3 eV 



and fa/C, > 2 x 10 7 GeV. 



(23) 



In the DFSZ model and grand unified models, C 7 ~ 0.75 so that m a $ 0.4 eV 



and f a > 1.5 x 10 7 GeV (Fig. pt). For models in which E/N = 2 and thus C 7 is 



very small, the bounds are significantly weaker. 

On the basis of their two-photon coupling alone, pseudoscalars can reach ther- 
mal equilibrium in the early universe. Their subsequent a — > 2j decays would 
contribute to the cosmic photon backgrounds [ p 



354 ], excluding a non-trivial 



m a -<? a7 -range (Fig. ]13|). Some of the pseudoscalars would end up in galaxies and 
clusters of galaxies. Their decay would produce an optical line feature that was 
not found 



387], |388j , p89| , leading to the "telescope" limits in Fig. |13j. For axions, 
the telescope limits exclude an approximate mass range 4-14 eV even for a small 



C 7 . 

Axions with a mass in the /ieV (10 -6 eV) range could be the dark matter of the 
universe (Sec. |6.3[ ). The Primakoff conversion in a microwave cavity placed in a 
strong magnetic field ( "haloscope" ) allows one to search for galactic dark-matter 
axions \ |355|1 . Two pilot experiments \ |390|, |391|1 and first results from a full-scale 
search 



392| ] already exclude a range of coupling strength shown in Fig. [13|. The 
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392, 



394] should cover the 



new generation of full-scale experiments 
dotted area in Fig. [T^, perhaps leading to the discovery of axion dark matter. 



6.2.2 Electrons 



7 + e - 
process e 



Pseudoscalars which couple to electrons are produced by the Compton process 
+ a [ 44, ^6], 49, |52| and by the electron bremsstrahlung 
(A, Z) -> (A, Z)+e~ + a [g|, g g |7(| [71], [74], [77|]. A standard 
solar model yields cLn cLxion luminosity of [ [48| L a — Oi ae 6.0 x 10 21 L Q where 
a ae is the axion electron "fine-structure constant" as defined after Eq. (|17|). The 



*ae ^ 2 X 10 



-23 



helioseismological constraint L a < 0.1 L of Sec. 2.3 implies a 
White-dwarf cooling gives [ ||, 39 a ae ^ 1.0 x 10~ 2B , while the most restrictive 
limit is from the delay of helium ignition in low-mass red- giants [ 52] in the spirit 
of Eq. (f 



a„ 



< 



0.5 x 10 



-26 



or 



9ae < 2.5 x 10™". (24) 

For m a > T ~ 10 keV this limit quickly degrades because the emission from a 
thermal plasma is suppressed. With Eq. (|l~7|) one finds for axions 



m, 



( C e < 0.003 eV and f a /C e > 2 x 10 9 GeV. 



(25) 



In KSVZ-type models C e = at tree level so that no interesting limit obtains. 
In the DFSZ model m a cos 2 (3 < 0.01 eV and f a / cos 2 (3 > 0.7 x 10 9 GeV. Since 
cos 2 j3 can be very small, there is no generic limit on m a . 
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6.2.3 Nucleons 

The axion-nucleon coupling strength is primarily constrained by the SN 1987A 



energy-loss argument [ jug fjjg, [m], g2g, g2§ [12J, [L2§ g2§, [L27|]. The main 
problem is to estimate the axion emission rate reliably. In the early papers it 
was based on a somewhat naive calculation of the bremsstrahlung process NN — > 
NNa, using quasi- free nucleons that interact perturbatively through a one-pion 
exchange potential. Assuming an equal axion coupling g a x to protons and neu- 
trons this treatment leads to the g a Ar-dependent shortening of the SN 1987A 
neutrino burst of Fig. ||. However, in a dense medium the bremsstrahlung pro- 
cess likely saturates, reducing the naive emission rate by as much as an order of 



magnitude [ 126]. With this correction, and assuming that the neutrino burst 
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Figure 14: Astrophysical and cosmological exclusion regions (hatched) for the 
axion mass m a , or equivalently the Peccei-Quinn scale f a . The globular-cluster 
limit depends on the axion-photon coupling; it was assumed that E/N = 8/3 
as in GUT models or the DFSZ model. The SN 1987A limits depend on the 
axion-nucleon couplings; the shown case corresponds to the KSVZ model and 
approximately to the DFSZ model. The dotted "inclusion regions" indicate where 
axions could plausibly be the cosmic dark matter. Most of the allowed range in 
the inflation scenario requires fine-tuned initial conditions. In the string scenario 
the plausible dark-matter range is somewhat controversial as indicated by the 
step in the low- mass end of the "inclusion bar." Also shown is the projected 
sensitivity range of the search experiments for galactic dark-matter axions. 
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was not shortened by more than half, one reads from Fig. || an excluded range 

3 x 1(T 10 < g aN < 3 x 1(T 7 . (26) 

With Eq. ([17]) this implies an exclusion range 

0.002 eV < m a C N < 2 eV and 3 x 10 6 GeV < f a /C N < 3 x 10 9 GeV. (27) 

For KSVZ axions the coupling to neutrons disappears while C p ~ —0.34. With a 
proton fraction of about 0.3 one estimates an effective Cjv ~ 0.2 so that [|4|, |126| 

0.01 eV < m a < 10 eV and 0.6 x 10 6 GeV < f a < 0.6 x 10 9 GeV (28) 

is excluded. 

In a detailed numerical study the values for C n and C p appropriate for the 
KSVZ model and for the DFSZ model with different choices of cos 2 (3 were im- 



plemented [ 127]. For KSVZ axions one finds a limit m a ^ 0.008 eV, while it 
varies between about 0.004 and 0.012 eV for DFSZ axions, depending on cos 2 f3. 
In view of the large overall uncertainties it is probably good enough to remember 



m a < 0.01 eV as a generic limit (Fig. 14) 



Axions on the "strong interaction side" of the exclusion range Eq. (gg) would 
have produced excess counts in the neutrino detectors by their absorption on 
oxygen if 1 x 10~ 6 < g a N i$ 1 x 1CT 3 [ |169||. For KSVZ axions this crudely 



translates into 20 eV < m a < 20 keV as an exclusion range (Fig. 14) 



6.2.4 Hadronic Axion Window 

This limit as well as the "trapping side" of the energy-loss argument have not 
been studied in as much detail because the relevant m a range is already excluded 
by the globular-cluster argument (Fig. [l4|) which, however, depends on the axion- 
photon interaction which would nearly vanish in models with E/N = 2. In this 
case a narrow gap of allowed axion masses in the neighborhood of 10 eV may 
exist between the two SN arguments ("hadronic axion window"). 

In this region one can derive interesting limits from globular-cluster stars where 
axions can be emitted by nuclear processes, causing a metallicity-dependent mod- 
ification of the core mass at helium ignition [ 51 1. It is intruiging that in this 



window axions could play a cosmological role as a hot dark matter component [ 
|395| 1 . Usually, of course, axions are a cold dark matter candidate. Moreover, in 
this window it may be possible to detect a 14.4 keV monochromatic solar axion 
line which is produced by transitions between the first excited and ground state 
of 57 Fe. In the laboratory one can then search for axion absorption which would 



give rise to x-rays as 57 Fe de-excites [ 396]. A recent pilot experiment did not have 



enough sensitivity to find axions [ 397], but a vastly improved detector is now in 



preparation in Tokyo (private communication by S. Moriyama and M. Minowa). 
6.3 Cosmological Limits 

The astrophysical axion mass limits are particularly interesting when juxtaposed 
with the cosmological ones which we thus briefly review. For f a > 10 8 GeV 
cosmic axions never reach thermal equilibrium in the early universe. They are 
produced by a nonthermal mechanism that is intimately intertwined with their 
Nambu-Goldstone nature and that implies that their contribution to the cosmic 
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density is proportional to /a' 175 and thus to m~ L175 . The requirement not to 
"overdose" the universe with axions thus leads to a lower mass limit. 

One must distinguish between two generic cosmological scenarios. If inflation 
occurred after the Peccei-Quinn symmetry breaking or if T re h ea t < fat t ne initial 
axion field takes on a constant value a\ = f a @i throughout the universe, where < 
0i < 7r is the initial "misalignment" of the QCD 6 parameter [ |39|, |39|, [4(1], |401 1 . 
If ©j ~ 1 one obtains a critical density in axions for m a ~ 1 /ueV, but since Gi is 
unknown there is no strict cosmological limit on m a . However, the possibility to 
fine-tune 0; is limited by inflation-induced quantum fluctuations which in turn 
lead to temperature fluctuations of the cosmic microwave background [ 402, 403 , 
404 , [405H - In a broad class of inflationary models one thus finds an upper limit 
to m a where axions could be the dark matter. According to the most recent 
discussion [ ^405f| it is about 10~ 3 eV (Fig. |i~4|) . 

If inflation did not occur at all or if it occurred before the Peccei-Quinn sym- 
metry breaking with T^he^ > f a , cosmic axion strings form by the Kibble mecha- 
nism [ 406 , 407 1- Their motion is damped primarily by axion emission rather than 
gravitational waves. After axions acquire a mass at the QCD phase transition 
they quickly become nonrelativistic and thus form a cold dark matter component. 
Unknown initial conditions no longer enter, but details of the string mechanism 
are sufficiently complicated to prevent an exact prediction of the axion density. 
On the basis of Battye and Shellard's treatment [ [408| , |409| and assuming that 
axions are the cold dark matter of the universe one finds a plausible mass range 
of m a = 6-2500 ^eV [|3|. Sikivie et al. [ 41C, 411, f412[| predict somewhat fewer 



axions, allowing for somewhat smaller masses if axions are the dark matter. 

Either way, the ongoing full-scale search experiments for galactic dark matter 
axions (Sec. |6.2.1| and Fig . |l~3| ) in Livermore (U.S. Axion Search [ |392| |) and in 
Kyoto (CARRAC K [ |393j , |394| | ) aim at a cosmologically well-motivated range of 
axion masses (Fig. |i~4| )^ 

7 LONG-RANGE FORCES 



7.1 Fifth Force 

New low-mass scalar or vector bosons would mediate long-range forces between 
macroscopic bodies. This is in contrast with pseudoscalars which couple to the 
spin and thus produce no long-range force between unpolarized bodies except 
for a residual force from two-boson exchange [ |413j , 414, |415| ]. In stars, a new 
long-range force has two different consequences. First, it modifies the effect of 
gravity. Second, it drains the star of energy, for the quanta of the new force are 
massless, or nearly so, and thus arise in thermal reactions. 



Thermal graviton emission is a case in point [ 416, 417, 415 , 419 , 420 |. However, 
the graviton luminosity is very small, about 10 -19 Lq for the Sun. Naturally, the 
coherent large-scale force is the most important aspect of gravity in stars! This 
conclusion carries over to new forces, notably a putative "fifth force." According 
to experiment, a fifth force has to be much weaker than gravity 
|42_, 

12 



424] so that possible modifications of stellar structure 
solar p-mode frequencies 



421, 422 



427, 



modifications of fundamental coupling constants near pulsars 
boson emission by the Hulse- Taylor binary pulsar 



425| , p6| or the 
are too small to be observable. Likewise, 

429] or scalar 



430 ] are negligible effects. 



However, there are no experimental fifth-force limits below about the centime- 
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ter scale, corresponding to boson masses exceeding about 10 -3 eV, where the 
most restrictive bounds arise from the energy loss of stars. The Yukawa cou- 
pling g$ (gv) of scalar (vector) bosons cf> to electrons has been constrained by the 
bremsstrahlung process e~ + 4 He — > 4 He + e~ + <fi which leads with Eq. (||) to 
limits of g s < 1.3 x 10~ 14 and gy < 0.9 x 10" 14 [ |, ||, |7]]. The Yukawa coupling 
to baryons has been constrained by the Compton process 7 + 4 He 
leading to g s < 4.3 x KT 11 and g v < 3.0 x 10~ n [ |, §6|, |7|]. 



l He + 



7.2 Leptonic and Baryonic Gauge Interactions 

It has been speculated that lepton and baryon number could play the role of 
gauge charges [ [43 1| , fi32|| . One consequence would be the existence of long-range 
leptonic and baryonic forces. The globular-cluster limits of the previous section 
translate into ej, < 1 x 10~ 14 and &b ~ 3 x 10~ n on the leptonic and baryonic 
gauge charges. Tests of the equivalence principle on solar-system scales constrain 
a composition-dependent fifth force, leading to something like eL,B ~ 10 -23 [ 
|433j ]. The cosmic background neutrinos would screen leptonic forces over large 
distances, but the solar-system limit on ei remains unaffected [ 1 13 [ . On the other 



hand, the SN 1987A neutrino burst would not suffer dispersion in the leptonic 
field of the galaxy because it is shielded by the cosmic background neutrinos [ 
|113f| . Leptonic forces contribute to the neutrino self-energy, modifying matter- 
induced neutrino oscillations in the Sun and supernovae [ 434, 435]. 



7.3 Time- Variation of Newton's Constant 

Astrophysics and cosmology are natural laboratories for testing all conceivable 
deviations from the standard theory of gravitatio n [ |436f[ . One hypothesis, which 
goes back to Dirac's large numbers hypothesis [ 437 , 438 1, holds that the value 
of Newton's constant Gn evolves in time. The present-day rate of change can be 
measured by a precision study of the orbits of celestial bodies. In the solar system 



data come from laser ranging of the Moon [ 436] and radar ranging of the planets, 
notably by the Viking landers on Mars [ 44C]. The increase of the length of day 
from 1663-1972 caused by tidal forces in the Earth-Moon system are consistent 
with a constant Gn 
Gn have been raised 

A weaker but less model- 



441], although some controversial claims for a decreasing 
436]. Beginning in 1974, very precise orbital data exist for 
the Hulse- Taylor binary pulsar PSR 1923+16 



4411 



dependent bound arises from the spin-down rate of the pulsar PSR 0655+64 [ 
445]. Finally, the long-time stability of galaxy clusters limits a decreasing Gn [ 
444 ] . The bounds from these methods are summarized in Table |2[ 

Very intruiging limits follow from the properties of the Sun. Paleontological 
evidence for its past luminosity, which scales approximately as Gj^Ai 5 , provides 
limits on previous values of Gn [ |445| ]. Solar models with a time- varying Gn 
were studied in the 1960s and 1970s \ E 



447, 448, 449, 45C, 451, 452], but 



truly interesting limits arose only recently from helioseismological observations 
A limit |Gn/G n | <2x 10" 12 — _1 



|453j , 454 , 455 ], A limit |Gn/Gn| ^ 2x 10~ i2 yr" 1 was derived by a comparison of 
the measured small-spacing p-mode frequency differences with those calculated 
from solar models with a time- varying Gn [ |455[|. The authors believe that the 



uncertainty is dominated by the observational errors while the prime systematic 
uncertainty is the exact solar age. A more conservative approach was used in 
deriving limits on a new solar energy loss mechanism (Sec. |2.3|) . The helioseismo- 
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Method G N /G N Authors Year Ref. 

[lO" 12 yr" 1 ] 





from 


to 










Laser ranging (Moon) 


-8 


8 


Williams et al. 


1996 




439 


Radar ranging (Mars) 


-12 


8 


Shapiro 


1990 




440 


Length of day 


-20 


20 


Morrison 


1973 




441 


Binary Pulsar 1913+16 





22 


Damour & Taylor 


1991 




442 


Spin-down PSR 0655+64 


-55 


55 


Goldman 


1990 




443 


Stability galaxy clusters 


-60 




Dearborn & Schramm 


1974 




444 


Helioseismology 


-2 


2 


Guenther et al. 


1998 




455 


Globular-cluster ages 


-35 


7 


Degl'Innocenti et al. 


1996 




461 


Pulsar masses 


-5 


4 


Thorsett 


1996 




462 



logical analysis of [ 455] provides the most restrictive limit on Gn/Gn, hopefully 



stimulating other groups to re-assess the bound independently. 

A large effect is expected on the oldest stars which "integrate" Gn(£) into the 



distant past. White-dwarf cooling is a case in point [ 45C, 457]. Under reasonable 



assumptions for the galactic age, the faint end of the luminosity function prefers 



a negative value for Gn/Gn around —10 to —30 x 10 -12 yr -1 [ 457]. Very re- 
cently this case has been re-examined [ |45§| 1 in greater detail. The observational 
uncertainty of the faint end of the luminosity function and the uncertainty of the 
galactic age preclude a clear limit on Gn/Gn- However, it is remarkable that 
even values as small as 10~ 14 yr -1 seem to make a noticeable difference for the 
cooling behavior of the oldest white dwarfs. 

Globular clusters are another important case because a different Gn in the past 
changes their apparent age based on the brightness of the main-sequence turn-off 



[ 459, 46C, 461]. A comparison with the expansion age of the universe brackets 
the allowed rate-of-change to the interval shown in Table |2[ 

A very sensitive limit arises from the observed masses of several old pulsars 
which measure the value of Gn at their time of formation in a SN explosion [ 
[462j ], The mass of a SN core at the time of collapse depends on its Chandrasekhar 
value which in turn scales as Gn 

The limits shown in Table [2] are difficult to compare on an equal footing as 
they involve vastly different ways of dealing with statistical and systematic un- 
certainties. However, it looks fair to conclude that |Gn/Gn| cannot exceed a few 
10 -12 yr -1 and that stars play an important role in this discourse. A similar 
bound arises from the cosmic expansion rate at the time of big-bang nucleosyn- 



thesis as measured by the primordial light-element abundances [ |463| . It implies 
that three minutes after the big bang Gn agreed with its present-day value to 
within a few tens of percent. A comparison of this limit with those of Table |2| 
requires a specific assumption about the functional dependence of G^(t). 

7.4 Equivalence Principle 



The general relativistic equivalence principle implies that the space-time trajec- 
tories of relativistic particles are independent of internal degrees of freedom such 
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as spin or flavor, and independent of the particle type (e.g. photon, neutrino). 
Several astronomical observations allow tests of this prediction. 

Limits on a gravitationally induced birefringence effect for photon propagation 
have been derived from the absence of depolarization of the Zeeman components 
of spectral lines emitted in magnetically active regions of the Sun [ 464]. Observa- 
tions of the light deflection by the Sun could soon become interesting [ 465]. The 
depolarization effect on distant radio galaxies already provides very restrictive 
limits 



468 



i 6q , |467| ] , as do pulsar observations 
One may also test for the equality of the Shapiro time delay between different 
particles which propagate through the same gravitational field. The absence of 
an anomalous shift between the SN 1987A photon and neutrino arrival times 
(Sec. |4.2| ) gave limits on violations of the equivalence principle [ 116, 117, 118 1 . 



|4.2| ) gave limits on violations of the equivalence principle 
The observation of a future galactic SN could provide independent arrival infor 
mation for v e 



lOg , [471] 



and v e and thus provide another such test 
A violation of the equivalence principle could manifest itself by a relative shift 
of the energies of different neutrino flavors in a gravitational field. For a given 
momentum p the matrix of energies in flavor space (relativistic limit) is E = 
p + M 2 /2p + 2p0(r)(l + F) where M 2 is the squared neutrino mass matrix, tp(r) 
the Newtonian gravitational potential, and F a matrix of dimensionless constants 
which parametrize the violation of the equivalence principle. F ^ can lead to 
neutrino oscillations in analogy to the standard vacuum oscillations which are 
caused by the matrix M 2 [ |47j| @74], |75|, [476], [477j, [47|, [47| 3 |48(], [481], |482 

14 



484]. Values for F$j in the general 10 —10 range could account for the 



solar neutrino problem. 



7.5 Photon Mass 

While it is usually taken for granted that photons are strictly massless, this 
theoretical expectation still needs to be tested experimentally. Some of the most 
restrictive constraints are related to the long-range nature of static electric or 
magnetic fields. The best laboratory limit of m 7 < 10 -14 eV derives from a test 
of Coulomb's law — see [ 485] for a review. 

In the astrophysical domain, the dispersion of the pulsed signal of radio pulsars 
is not a very sensitive diagnostic as the interstellar medium mimics a photon mass 
corresponding to a plasma frequency of order 10 -11 eV. The spatial variation of 
magnetic fields of celestial bodies is far more sensitive. Jupiter's magnetic field 

" 15 eV [ ||| while the Earth's 

I0~ 15 eV [|87]]. 



as measured by Pioneer-10 yields m 7 < 0.6 x 10 
field gives m 7 <0.8x 



If the photon has a mass, Ami is an observable quantity where A is the 



vector potential corresponding to known magnetic fields, 
experiment discloses Am 2 < 0.8 x 10~ 22 T m eV 2 [ || 
field implies A~2x 10 9 T m (Tesla-meter) so that m 7 <2x 
cluster-level field corresponds to A ~ 10 12 T m, providing va. 



A recent laboratory 
The galactic magnetic 
10~ 16 eV while a 
10~ 17 eV. 



< 

7 r^i 



Even more restrictive limits obtain from astrophysical objects in which mag- 
netic fields, and hence the Maxwellian form of electrodynamics, play a key role 
at maintaining equilibrium or creating long-lived stable structures [ |489| ]. The 
most restrictive case is based on an argument about the magneto-gravitational 
equilibrium of the gas in the Small Magellanic Cloud. The argument requires 
that the range of the interaction exceeds the characteristic field scale of about 
3 kpc [ 49C] . This resulting limit m 7 < 10~ 27 eV, if correct, is surprisingly close to 



Particle Physics from Stars 



37 



10 33 eV where the photon Compton wavelength would exceed the radius of the 
observable universe and thus would cease to have any observable consequences. 



7.6 Multibody Neutrino Exchange 

Two-neutrino exchange between fermions gives rise to a long-range force. A 
neutrino may also pass around several fermions, so to speak, producing a much 
smaller potential. In a thought-provoking paper it was claimed that this multi- 
body neutrino exchange could be a huge effect in neutron stars, essentially be- 
cause combinatorial factors among many neutrons win out against the small- 
ness of the potential [ 491]. To stabilize neutron stars, it was claimed, the 
long-range nature of neutrino exchange had to be suppressed by a nonvanish- 
ing mass exceeding about 0.4 eV for all flavors. In an interesting series of pa- 
pers it was shown, however, that a proper resummation of a seemingly divergent 
series of terms leads to a well-behaved and small "neutron-star self-energy" [ 

limit. 



493, 494, 495, 496, 497], invalidating the claim of a lower neutrino mass 



8 CONCLUSION 

Stellar-evolution theory together with astronomical observations, the SN 1987A 
neutrino burst, and certain x- and 7-ray observations provide a number of well- 
developed arguments to constrain the properties of low-mass particles. The most 
successful examples are globular-cluster stars where the "energy-loss argument" 
was condensed into the simple criteria of Eqs. (||) and @ and SN 1987A where it 
was summarized by Eq. dH). New particle-physics conjectures must first pass these 
and other simple astrophysical standard tests before being taken too seriously. 

A showcase example for the interplay between astrophysical limits with labora- 
tory experiments and cosmological arguments is provided by the axion hypothesis. 
The laboratory and astrophysical limits push the Peccei-Quinn scale to such high 
values that it appears almost inevitable that axions, if they exist at all, play an 
important role as a cold dark matter component. This makes the direct search 
for galactic axion dark matter a well- motivated effort. Other important standard 
limits pertain to neutrino electromagnetic form factors — laboratory experiments 
will have a difficult time catching up. 

The globular-cluster limit was based on relatively old observational data. A 
plot like Fig. || could be made more significant with dedicated CCD observations 
of globular clusters and improved theoretical interpretations. Assuming that such 
an effort produces internally consistent results, the statistical significance would 
improve, but I would not expect a vast gain for, say, the neutrino magnetic- 
moment limit as there always remain irreducible systematic uncertainties. 

Shockingly, SN 1987A as a particle-physics laboratory is based on no more than 
two dozen measured neutrinos. The observation of a future galactic SN with a 



large detector like Superkamiokande or a future observatory such as OMNIS [ 247] 
would provide a high-statistics neutrino light curve and thus a sound empirical 
basis for SN theory in general and for particle-physics interests in particular. 
Alas, galactic supernovae happen only once every few decades, perhaps only once 
per century. Thus, while the neutrinos from the next galactic SN surely are on 
their way, it could be a long wait until they arrive. 

Most of the theoretical background relevant to this field could not be touched 
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upon in this brief overview. The physics of weakly coupled particles in stars is 
a nice playing field for "particle physics in media" which involves field theory 
at finite temperature and density (FTD), many-body effects, particle dispersion 
and reactions in magnetic fields and media, oscillations of trapped neutrinos, and 
so forth. It is naturally in the context of SN theory where such issues are of 
particular interest, but even the plasmon decay 7 — > vv in normal stars or the 
MSW effect in the Sun are interesting cases. Particle physics in media and its 
astrophysical and cosmological applications is a fascinating topic in its own right 
which well deserves a dedicated review. 

Much more information of particle-physics interest may be written in the sky 
than has been deciphered as yet. Other objects or phenomena should be consid- 
ered, perhaps other kinds of conventional stars, perhaps more exotic phenomena 
such as 7-ray bursts. The particle-physics lessons to be learned from them are 
left to be reviewed in a future report! 
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